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Results of the idealized mode-coupling theory for the structural relaxation in suspensions of 
hard-sphere colloidal particles are presented and discussed with regard to recent light scattering 
experiments. The structural relaxation becomes non-diffusive for long times, contrary to the ex- 
pectation based on the de Gennes narrowing concept. A semi-quantitative connection of the wave 
vector dependences of the relaxation times and amplitudes of the final a-relaxation explains the 
approximate scaling observed by Segre and Pusey [Phys. Rev. Lett. 77, 771 (1996)]. Asymptotic 
expansions lead to a qualitative understanding of density dependences in generalized Stokes-Einstein 
relations. This relation is also generalized to non-zero frequencies thereby yielding support for a 
reasoning by Mason and Weitz [Phys. Rev. Lett 74, 1250 (1995)]. The dynamics transient to the 
structural relaxation is discussed with models incorporating short-time diffusion and hydrodynamic 
interactions for short times. 
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I. INTRODUCTION 



The dynamics of suspensions of colloidal particles has 
been the topic of active research for many years |l|; ^ . 
Whereas the motion of isolated Brownian particles has 
been well understood for long, less is known about the 
dynamics of concentrated suspensions. Direct particle 
interactions and solvent mediated hydrodynamic interac- 
tions (HI) are important if the colloidal volume packing 
fraction increases above a few percent |l] . Experimental 
studies mainly employing dynamic light scattering (DLS) 
have provided a wealth of information on dense systems 
and are also the stimulus for the theoretical work pre- 
sented in this contribution. 

The pioneering study by van Megen and coworkers 
of the liquid to glass transition in hard-sphere like col- 
loidal dispersions has provided detailed data on the den- 
sity fluctuations at this dynamic, ergodic-to-nonergodic 
transition ||^-[ictl. Besides their intrinsic interest, these 
experiments also made possible quantitative tests ]^pT| 
of predictions from the idealized mode-coupling theory 
(MCT) ||l^,|T^. Agreement of experiment and theory 
within errorbars of ca. 15 % has been reported. This 
comparison, which up to now has tested leading order 
asymptotic predictions and has thus restricted the valid- 
ity of the theoretical results to small separations from the 
critical density, provides support for the glass transition 
scenario as described by MCT, which has also been stud- 
ied for colloidal micronetwork spheres charged 
colloidal systems and colloidal emulsions |l^]. Re- 
cent DLS experiments by Segre, Pusey and coworkers 
study hard sphere like systems at lower colloidal densi- 
ties and report unexpected and seemingly unrelated scal- 
ing properties of the dynamic scattering functions 
pO| . Thus the question arises for which density range 
below the glass transition the MCT describes the dom- 



inant physical mechanism observed in the dynamics of 
concentrated colloidal fluids and whether the reported 
scalings can be explained by MCT. Studies of a gener- 
alized Stokes-Einstein relation |19; 21| and optical mea- 



surements by Mason and Weitz p2[ further raise the 
question about the connection of the viscoelastic moduli 
p3[ to the collective and self particle motion at rather 
high densities which can also be considered using MCT 

The MCT was developed starting from theories of the 
dynamics of simple liquids upon the realization that in 
that subsystem of the equations of motion which aimed 
at describing the structural relaxation there exists a bi- 
furcation separating ergodic from non-ergodic motion 
[^,|2^. The physical mechanisms held responsible have 
been called "cage effect" and "back flow" phenomenon 
[|l3|,|2^. This transition was suggested as origin of the 
slowing down and of the anomalies of the dynamics at the 
glass transition. The idealized MCT studies the struc- 
tural relaxation neglecting all other, possibly present, 
long-time dynamical effects ||l2 , 
MCT discusses long-time ergodi 
tions [E7|,|30|-|2~ 
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The extended 
icity restoring correc- 
The bifurcation at critical values of 
the thermodynamic parameters like the colloidal packing 
fraction, 0, introduces a small (separation) parameter, 
e — {(j) — (j)c)/4>c, and the possibility of asymptotic expan- 
sions in e; see Refs. [^^,|3j] for references and detailed 
results. 

Two asymptotic scaling law regions can be shown. In 
the first, for intermediate times, the feedback mechanism 
of caging of particles, causes an ergodic-to-nonergodic 
bifurcation, which is characterized by universal power 
law decays. During the second, for longer times, the col- 
lective rearrangements of the cages requires cooperative 
dynamics, like the build up of back flow patterns first 
discussed for liquid helium p3-p^. The strongly corre- 
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lated dynamics manifests itself in a coupling of the time 
scales for this final process of the structural relaxation. 
As it describes, in the liquid, the decay of the incipient 
frozen glassy structure, it is not surprising, that its MCT 
description requires detailed information about the equi- 
librium structure. 

As the leading-in-e asymptotic results exhibit numer- 
ous nontrivial universal features, experimental tests of 
MCT mainly address these and thus a number of correc- 
tions need to be considered: (z) The structural relaxation 
itself leads to corrections of higher order in e which limit 
the range of validity of the leading asymptotics. (m) The 
short-time or microscopic dynamics affects the transient 
to the structural dynamics and needs to be considered if 
no clear separation of time scales is possible, (iii) Long- 
time relaxational mechanisms may be present, which by- 
pass the structural relaxation, and lead to faster decay. 
The third correction appears to be absent in colloidal sus- 
pensions at the densities of interest for the present study 
and thus shall be neglected in the following It 
is interesting to mention, though, that in colloidal emul- 
sions droplet shape fluctuations cause long-time relax- 
ation and can be explained within the extended MCT 
|p7| . Theoretical understanding of the first correction 
effect in lowest relevant order in e has been achieved 
recently and is the basis of the present consid- 

erations. Lacking a deeper understanding of the micro- 
scopic transport effects of colloidal suspensions (point {ii) 
above) a qualitative study shall be undertaken incorpo- 
rating short-time effects with the most simple approx- 
imations compatible with the MCT description of the 
structural relaxation. Thus the limit of the predominance 
of the structural relaxation is estimated from numerical 
solutions of the MCT equations using simple models for 
the microscopic transient dynamics of colloidal suspen- 
sions. Brownian short-time diffusion with and without 
hydrodynamic interactions is considered 

The paper is arranged as follows. Section II summa- 
rizes the equations of motion of the idealized MCT. Sec- 
tion III presents and discusses our results, focusing first 
on the aspects purely structural-relaxational and then 
on the microscopic transient effects. The mentioned ex- 
perimental findings are addressed in section III.C. Short 
conclusions end the paper. 



II. EQUATIONS OF MOTION 



tral quantity, the (normalized) intermediate scattering 
function, ^q{t) — ■^{g*{t)gq{0))/Sq, measures the time 
and wave vector dependence of the collective density fluc- 
tuations, gq{t), around the average homogeneous density 
which for hard-sphere particles of diameter d is converted 
to the volume packing fraction. The normalization is pro- 
vided by the static structure factor, Sq — {\gq\'^)/N 

The equations of motion of the idealized MCT from 
which asymptotic analysis extracts the physically rele- 
vant long-time dynamics can be summarized as 

^q{t) - mq{t) + i dt' rUqit - t' ) $,(t') = , (1) 
"^9 W = Ek+p=q ^(9; k,p) $fc(i) %{t) , (2) 

where the coupling constants or vertices, V ^ are uniquely 
specified by the static structure factor, Sq] see Refs. 
P6|,|39[| for explicit formulae. Equation (p|) approximates 
the autocorrelation function of the fluctuating forces by 
considering a force to arise between two density fluctu- 
ations interacting via an effective potential. Then the 
correlation function of the four density fluctuations is 
approximated by the squared density correlators and the 
effective potential enters the vertices. 

In Eqs. (0J|) the structural dynamics results from the 
equilibrium state of the fluid as captured in Sq which is 
the only input. Thus the long-time structural relaxation 
of a dense suspension of interacting Brownian particles is 
predicted to be identical to the one of an atomic system if 
the interaction potentials of both systems coincide. Here 
we will consider hard sphere potentials only. By ansatz, 
MCT neglects the possibility of an ordered, crystalline 
state and thus in experimental comparisons crystalliza- 
tion has to be prevented. Then the structure factor Sq of 
a liquid of hard spheres is known to be a regular function 
smoothly varying with packing fraction |38| and conse- 
quently the vertices in Eq. (0) are regular functions of 
the single (for a liquid of hard spheres) thermodynamic 
state parameter (j). 

Brownian particles diffuse on (appropriately chosen) 
short distances and thus for short times [Q. Incorporat- 
ing this into MCT leads to the following simple model of 
colloidal suspensions close to the glass transition ^ : 



$,(t)-m,(<) + - 
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dt' mq{t-t')<^q{t') = 



(3) 



The idealized MCT of the liquid to glass transition 
leads to a closed set of non-linear equations of motion 
for the density fluctuations jlj, ^|2^. Other dynamical 
variables are connected to them via the Zwanzig-Mori 
formalism and mode-coupling approximations. The the- 
ory aims at a description of the structural relaxation as 
it emerges from the (microscopic) short-time dynamics 
and slows down owing to the increasing density and thus 
increasing importance of particle interactions. The cen- 



This equation replaces Eq. (|l|) and extends it to short 
times where the initial condition runs $q(t) = l — q^D'^ t. 

Two approximations for the short-time diffusion coeffi- 
cient, Z?^' , are widely used in theoretical work on colloidal 
dispersions and differ in the treatment of the solvent ef- 
fects. In the most simple model of Brownian diffusion 
the interaction of the solvent with the colloidal particles 
is modeled with a single friction coefhcient, Co This 

leads to Dq'^^^ = Dq/ Sq, where Dq is given by Einstein's 
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law, Do = ^bT/^o, and the denominator arises from par- 
ticle interactions as was first argued by de Gennes in a re- 
lated context . This approximation is not satisfactory 
except for very low packing fractions, because the solvent 
also leads to long-ranged and quasi-instantaneous inter- 
actions of the colloidal particles, called hydrodynamic in- 
teractions (HI). Whereas the HI do not affect the equi- 
librium statistics and thus Sq of the colloidal system, 
their effects on short time scales cannot be neglected in 
general and are captured in a wave vector dependent am- 



plitude, Hq 



n 



s.(HI) 



DoHq/Sq. Progress on a de- 



tailed theory for Hq has proven very difficult but has 
culminated in accurate results for it up to intermediate 
packing fractions ^Sj. At the considered high pack- 
ing fractions however, Dq''™^ can only be estimated from 
experiments or simulations at present |18, 4^. Note, that 
our approach to incorporate HI into the MCT equations 
of motion only via D^' differs from the one developed 
by Nagele and others jQ, which aims at describing the 
dynamics at lower packing fractions, and which would 
affect the structural long-time dynamics. Our approach 
also differs from the work of Cohen et al. who incorpo- 
rate aspects of the cage effect into an effective short-time 
Dg 0,|7). The role of the HI here also differs from a 
recent theory of Tokuyama et al. ||48|-|50[| who consider 
the HI in noncquilibrated colloidal suspensions. 

Equations (^) and (||) have been solved repeatedly 
without HI and with different approximations for the 
structure factor of the hard-sphere fluid p^-p^, Q ; for 
the details of the numerical calculations see the quoted 
references. Various aspects of the known solutions will 
be connected to recent experimental observations in this 
contribution and new solutions taking HI into account 
via Dq^^^^ will be presented, which are an extension of 
the calculations in Refs. 

Figure |l] shows the short-time diffusion coefficients en- 
tering the numerical calculations discussed in the fol- 
lowing. The short time diffusion coefficient without HI, 
follows immediately from the hard sphere struc- 
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ture factor, where the Percus-Yevick approximation is 



used 



The Sq shown also enter the vertices in Eq. 



The short-time diffusivity with HI, U, 



■(HI) 



is cho- 



sen as shown in Fig. |l|. It is aimed at a discussion of the 
dynamics with HI transient to the structural relaxation, 
and thus, for the high densities considered, a rough ap- 
proximation modeled from the experiments in is 

used. Outside the window 0.77 < q/qp < 2.4, Dq is 

assumed constant for simplicity. The values of Dq^™^ for 
q ^ 0, q = qp and q ^ oo are adjusted to 0.2:1.74:1.0 
mimicking the measured ratios ||l8|,Q|. Within the men- 
tioned wave vector window, the experimental data are 
modeled by: D° / Dq^™^ = — x(g/gp)^/ In/^^^^, where 
is the MCT critical non-ergodicity parameter, and 
X = 0.29 leads to a continous matching. "De Gennes 
narrowing" is present in D^' in both approximations, as 
its inverse varies in phase with the structure factor 5*^ 



and an appreciable slowing down of the short-time dy- 
namics for wave vectors q around the principal peak at 
q = qp results. 

Some representative numerical results for the collec- 
tive density correlators, ^q{t), obtained as specified in 
are exhibited in Fig. ||, where also the shown 
wave vectors are indicated. The correlators of Fig. ||(a) 
correspond to a density rather close to the critical liquid- 
to-glass bifurcation point, ipc — 0.516, of this model. The 
reduced distance equals e = {(/)—(f)c)/4'c — — 10~"/^, with 
n = 9, where, as in the following, in order to simplify 
comparison with Refs. ||3^, the packing fractions will 
be reported by stating the number n. Decreasing the 
packing fraction to e = —0.1 corresponding to n — 3, 
results in the intermediate scattering functions of Fig. 
^(b). Results with and without HI as shown in Fig. g 
for various packing fractions will be discussed in the fol- 
lowing. Only the dynamics in the colloidal liquid phase 
is shown, where the correlators decay to zero during the 
final relaxation process, because the mentioned experi- 
mental studies focus on this final decay; for MCT results 
on hard sphere glasses at higher packing fractions see 
Refs. [^,^, It is also interesting to note that short 
range attractions can increase the colloidal glass packing 
fraction appreciably ]5l|,[52|. 



III. RESULTS AND DISCUSSION 

A. Leading asymptotic scaling laws 

In lowest order in the separation parameter e, the 
MCT predicts the existence of two divergent time scales 
with two different scaling laws describing the dynamics 
in expanding windows in time or frequency; see the Refs. 
||l2|,|2^, ^^,|5^ for detailed derivations and reviews of 
these results. 

In the first or /^-scaling law window, a factorization 
property allows to separate the sensitive and rather uni- 
versal dependences of the dynamics on the separation 
parameter and on time from the system specific depen- 
dences like spatial variation. 

%it) - fq + hq G{t, e) for \%{t) - « 1 . (4) 

The /^-correlator is given by a homogeneous function, 
G cx \/\£\ g^(t/t^), specified by one system specific pa- 
rameter A, which can be calculated for simple liquids 
from Sq and determines all exponents of MCT |^ . The 
first divergent scaling time ~ to lel^^^/^"-*, lies in the 
center of the window of validity of Eq. (Q) and, below 
the critical density, can be taken from the root of G: 
^q{t — tg) — f^. The one parameter to, the crossover 
or matching time, remains as only remnant of the short- 
time or transient motion and can only be obtained from 
matching the asymptotic results to the full dynamics in- 
cluding some short-time model. In Fig. |^(a) one notices 
that to differs by a factor 1.2 for the two models of the 
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transient. A shift of the curves with HI relative to the 
ones without HI collapses both sets of curves for times 
t^O.lcP / Dq. The great simplification of the dynamics 
provided by the factorization in Eq. may be inter- 
preted as resulting from a localization transition close to 
which density fluctuations relax via local rearrangements 
and not via mass-transport over larger distances. If the 



spatial variation of and hq is studied in detail 



the 



localization may be traced back to the "cage-effect" , that 
particles are surrounded by next-neighbor shells whose 
ability to cage the central particles depends on the fluc- 
tuations of the local structure and thus, in a cooperative 
manner, on the dynamics of the caged particles them- 
selves. 

In the second scaling law region, another set of diver- 
gent time scales, called a-relaxation times Tq, appears 
and a-master curves describe the final relaxation of the 
density correlators from fq to zero during that time win- 
dow fB3|- 



percooled atomic liquids [Q. Although the collective 
density fluctuations of the colloidal Brownian particles 
are diffusive on short time scales due to random colli- 
sions with solvent molecules, during the structural relax- 
ation only stress fields arising from colloid-colloid par- 
ticle interactions survive out to long times. Thus large 
distance density fluctuations decay by local particle rear- 
rangements. The strong slowing down of *&g(i) on length 
scales of the order of the average next -neighbor distance 
indicates local and cooperative particle rearrangements 
and is reminiscent of the back flow phenomenon familiar 
from simple liquids ||35|-|37|] . 

The coupling of the wave vector modes in Eqs. (|^J^) 
explains the qualitative trend that the correlators with 
larger g-process amplitudes, J^, relax slower, i.e. have 
a larger Tq 
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|33[ |. Intriguingly, for a hard sphere 
liquid the wave vector dependence of the dimension-less 
time scale Tq'' — —r^/(cPhifq), is rather close to the 
one of the actual a-time, Tq, at least for intermediate 
wave vectors. Note that the comparison shown in Fig. 



<^q{t) -> <^q{t/Q for e~^Q- , t/t'^= fixed . (5) | must be taken with a grain of salt, as the definition 



This superposition principle states that the final relax- 
ation processes (asymptotically) depend on the distance 
to the critical point only via the relaxation times Tq, 
which moreover are coupled, Tq = fqt'^, and diverge upon 
approaching (/)c: t'^ = to|e|~''' with 7 > (l/2a). The equa- 
tions, which the obey, are obtained in a special limit 
from Eqs. (|l|,^ and are consequently independent of the 
microscopic short-time dynamics The resulting 

two-step relaxation scenario of the idealized MCT has 
been fully worked out for some simple liquids; see Refs. 
1^,^, ^ for calculations of the exponents and master 
functions of the two scaling regions for a hard sphere 
liquid. 

The quantities of most immediate interest to experi- 
mental observations of the final or a-relaxation process 
are the relaxation times |5^. Figure ^ presents results 
for the asymptotic Tq from the model specified above 
and compares them to previous calculations using a dif- 
ferent approximation (Verlet-Weis form) for the static 
structure factors of a hard sphere liquid . Very small 
differences in the Tq result from the two approximations 
to Sq. A discussion of short-time sum rules for colloidal 
suspensions as done in the spirit of de Gennes ^] leads 
to the prediction of (short-time) relaxation times obey- 
ing Tq — (l/q^Dq). Such a behavior for the Brown- 
ian model, scaled to match Tq for q — qp, is also indi- 
cated in Fig- 0- The MCT a-relaxation times obtained 
from Eqs. (^JJ), where the transient does not enter, and 
the results from the short-time sum rules qualitatively 
are similar for not-too-small wave vectors because both 
vary in phase with the structure factor. Their differ- 
ent physical origins, however, clearly show up for small 
wave vectors where the short-time relaxation times, Tq , 
become diffusive, whereas the MCT a-relaxation times 
Tq become wave vector independent as first anticipated 
in Mountain's description of Brillouin scattering in su- 



of Tq is not unique because of the stretching, i.e. non- 
exponentiality, of the a-process in MCT. Nevertheless, 
this semi-quantitative connection of the a-process am- 
plitude to its time scale, suggests a possible (partial) col- 
lapse of the ^q{t) for different q at the same packing 
fraction onto a common curve given by: 



Ar^ ft) 

$,(i) = cxp{-— :^}, 



6d2 



where ri^-* — 



(6) 



Conceptually, Ar'^{t) should be connected to the mean- 
squared displacement of a colloid particle, to be denoted 
by Sr'^{t). From the definition of Tq \ Eq. (||), and the 
factorization property, Eq. (^) , immediately follows that 
very close to the critical packing fractions all rescaled 
curves intersect at the /3-scaling time t^. The connec- 
tion of fq via Tq'' to the a-relaxation time, see Fig. ||, 
then shepherds the correlators to stay close during the 
final relaxation step, too. Figure ^ shows representative 
scaled correlators for two packing fractions, where the 
used wave vectors are marked in the inset. The correla- 
tors are drawn for $q(t) > 0.05 in order to prevent over- 
crowding the figure. The g-dependent stretching of the 
correlators causes a noticeable spreading of the rescaled 
correlators for long times. Considering Fig. ^, one also 
does not expect a data collapse for wave vectors outside 
the shown q-range. Moreover, this scaling explicitly vio- 
lates the short time behavior of the intermediate scatter- 
ing functions which e. g. become diffusive for small wave 
vectors invalidating Eg. (|^). This explains the spread 
of the curves in Fig. ^ at short times. In Fig. ^ also 
the mean-squared displacement from Ref . 1 34 1 shifted as 
suggested by Eq. (||) is shown and lies within the clatter 
of the curves. 

This ansatz, together with the known Gaussian ap- 
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proximation to the self intermediate scattering function 
[0,^, ^ , gives a most simplistic description of the coher- 
ent and incoherent density correlators of the MCT. Nev- 
ertheless, the only point where Eq. (|^) asymptotically 
rigorously collapses all correlators is at ^qite) — fg be- 
cause of the factorization property, Eq. (Q). Already in a 
vicinity of this point a spread of the curves exists because 
of hq hMSofq/idd'^Tq''), which would follow from Eq. 
and the known /3-expansion, Sr^{t)=r^^ — husoGit) 




B. Corrections 

The discussion up to now has used the asymptotic for- 
mulae to lowest orders in the separation parameter and 
thus might restrict the discussed phenomena to close 
neighborhoods of the critical packing fraction 0c- The 
leading corrections in e to the asymptotic scaling laws 
of Sect. Ill A have recently been discussed in detail for 
the present model [^,|3j], and in some cases allow to ex- 
tend the range of validity of the asymptotic expansions 
appreciably. 

The corrections to the /3-scaling law, Eq. (|^), for the 
dynamics close to are of the form: 

- + hq ( G{t) + H{t) + Kq G\t) +eKq), (7) 

where the Kq and Kq are wave vector dependent con- 
stants which follow from asymptotic solutions to Eqs. 
(|l|,^. See for the definitions and for the correc- 

tion function H{t), which is of order 0{e). The range 
of validity of the /3-scaling law, Eq. (^, is thus found 
to be of order ^/e, and to differ for different wave vec- 
tors or observables. The /3-region description of Eq. (j^) 
extends the range of usefulness of the MCT asymptotic 
expansion around the critical non-ergodicity plateau ap- 
preciably as can be seen in Refs. p3|,|3^, and provides 
detailed few parameter formulae for the density correla- 
tors which have already found use in the data analysis 
of computer simulation studies [p8|-^6l[. For the curves 
without HI of Fig. ^, Eq. (0) describes the correlators in 
the /3-window on a 10% error level starting from the time 
scale t^Q.l(P /Dq which was estimated in Sect. 



Ill A to 



be the range of domination of the structural relaxation. 
Thus Eq. (Q) extends the asymptotic expansions of the 
structural relaxation almost to the microscopic dynamics 
in this model. 

The Qf-process has been the focus of the recent DLS 
scattering studies |Q and, for wave vectors around the 
peak, describes the main portion of the decay of ^q{t). 
The range of validity of the a-process-superposition 
principle, Eq. (|^), is appreciably larger than the one 
of the /3-scaling law: 

^ ^qit/Q + e^qit/Q for e ^ - . (8) 

The corrections in Eq. (||) are only of linear order in e. 
Although the complete form of '^q{t) is not known yet. 



its variation for short rescaled times can be deduced and 
can be argued to give the dominant correction to Eq. (j^) 
for not too large e: ^q{i ^ 0) = —hq Bi t^^, where the 
coefficient Bi is of order unity. As this term can grow 
without bounds, the dominant aspect of the leading cor- 
rections is to cause the correlators (t) to rise above the 
Qf-master curves for times shorter than the a-relaxation 
time. In this time window, the a-master curves follow 
von Schweidler's law, <&g(<) ^ fq = —hqBr (again with 
B — 0(1)) 1^,^. As the g-dependence of the time 
scales of the a-correlators $q(t) can be estimated from 

^S)=m - (V^r^)'), with rf^) = {Brjhq)^/»t', 
the short time corrections can be rewritten in the 



time window t 



< t < i': 



(vS)' 



r(vS)N,-6- 



(9) 



Thus it is apparent that deviations from the asymptotic 
a~process scaling law, Eq. (^, are stronger for corre- 
lators with a shorter a-relaxation time or smaller a- 
process amplitude; the second connection arising because 
of the relation between Tq and see Fig. |[ 

If at larger separations from the critical density the 
Qf-relaxation times are determined from the correlators 
^q{t), then the corrections to the a-scaling law, Eq. (||), 
may differently affect Tq{<j>). This is caused by the in- 
herent stretching in the a-master curves 1 32 and by the 



time variation of the corrections ^g(t). The dominant 
short time variation of ^q(t) leading to Eq. (^) will af- 
fect Tq {(/)) if a definition of the relaxation times is used, 
which stresses the initial decay during the a-process. A 
possible definition of Tq{(j)) exhibiting this effect is given 
by: ^o{t = Tq) — \fq- Some results are indicated in 
Fig. 1^, where also the a-process amplitudes, are 
shown in the inset. As Fig. || shows, this definition of Tq 
asymptotically gives almost identical g-dependences as 
obtained from the ^q{i) 

Because of the rather large range of validity of the a- 
scaling law, Eq. (|^), for the intermediate scattering func- 
tions as explained by Eq. (|^) , only very small deviations 
of rq(0) from the asymptotic wave vector dependence are 
seen in Tq{—e)'' for n > 3. As expected from Eq. (||), 
for 71 = 3, which lies close to the limit of applicability 
of a-scaling, the relaxation times are relatively longer 
and the largest (relative) deviations appear for correla- 
tors with small or Tq. At this separation from the crit- 
ical density, (j) = 0.9 0c, already some differences for the 
two models of the short time diffusion, with and without 
HI, are noticeable in Fig. As shown in the inset, the 
differences can almost completely be incorporated into a 
packing fraction dependent shift of the matching time to- 
ll the time scales are normalized to unity for q = Qp, then 
collapse can be achieved of the at n = 3 except for the 



5 



smallest wave vectors. Note that some finer aspects of 
the figure depend on the special choice how to measure 
Tq. For example, the correlators without HI at n = 3 
and for q = Qp and q = 0.9Aqp (just below it) actually 
almost overlap and the apparent differences in Tq arise 
solely from the /^-values entering the used definition. 

For even larger separations from the critical density, 
n = 2 and n = 1 in Fig. ||, clear differences of the 
long-time scales with and without HI appear and can 
obviously not be explained by structural relaxation, Eqs. 
(|[||), alone. The diffusive particle motion on short time 
scales causes the correlators for small wave vectors to 
decay slower relative to the non-diffusive a-process. 

Often a diffusive behavior is assumed also for the struc- 
tural relaxation and the relaxation times are converted 
to diffusion coefficients via 1/Dq ~ q^Tq. Fig. ^ shows 
so calculated Dq normalized aX q — qp in order to elim- 
inate the drift of qp with packing fraction; see Fig. |^. 
Almost no deviations from the asymptotic variation as 
follows from the a-scaling law, Eq. (^, can be be seen 
for n > 3. Note that the non-diffusive character of the 
structural relaxation is hidden in this representation. For 
larger separations and thus smaller packing fractions, a 
trend of the long-time diffusion coefficients with hydro- 
dynamic interactions (HI) to approach the shape of the 
short time ones, D^q'^™^ , can be recognized. 

Considering Figs. || and ^ one needs to keep in mind, 
however, that differing methods to determine the final 
relaxation times or the long-time diffusion coefficients 
would lead to somewhat different g-dependences because 
they would weigh the stretching of the a-process ^q{i) 
and the leading corrections ^'g(i) differently. The defini- 
tions chosen here allow to explain the wave vector and 
packing fraction dependences in Tq and Dq from known 
aspects, Eqs. (H,^), of the asymptotic expansions. 



C. Visco— elastic properties 

The time or frequency dependent shear modulus, G^, 
of colloidal suspensions can be defined as an autocor- 
relation function of elements of the stress tensor and 
splits into three contributions |6|, ^ . Whereas the first 
arises from the direct potential interactions of the parti- 
cles and is familiar from simple atomic liquids, the lat- 
ter two contain effects of the HI and are peculiar for 
colloidal particles immersed in a solvent. Only for the 
first potential part there exist MCT expressions which 
are applicable close to the glass transition at (j)c 
^,|3|]; however see Ref. |Q for lower densities. Sim- 
ilarly as for mq(t) from Eq. (^), G,,(t) is given by a 
quadratic mode-coupling functional in the $g(t). Con- 
sistent with the neglect of the HI contributions to (t) , 
solutions for the $g(i) are used which are calculated 
without HI, i. e. with the short-time diffusion coeffi- 
cients Dq^'^K Figure]^ shows the frequency dependent 
storage and loss shear moduli for a number of densities 



[Q. As we consider the part of Grj arising from po- 
tential colloidal interactions only, and thus cannot ad- 
dress the importance of HI at higher frequencies, only 
results in the frequency window of structural relaxation 
are shown. For low frequencies, the viscosity rj can be 
obtained via GJJ((jJ — > 0) — > ujii] — ?/oo), where rjoo is 
the high-frequency shear viscosity which is caused by in- 
stantaneous solvent interactions |^, We use the 
approximation ijoo = ksT/ (SndDo). A plateau region 
in G^(w) corresponds to the /3-scaling window, Eq. (||), 
and indicates elastic behavior of the colloidal suspension 
on intermediate time scales. In the non-ergodic states 
above (j}c, the colloidal system would be characterized by 
a finite elastic shear modulus, Gjj> C^, where the value 
at the glass transition follows from the fq. The appro- 
priate Fourier transforms of the /3-correlator describe the 
dynamics around this elastic plateau in GJ^, and in the 
minimum region of G"(a;) between the transient high fre- 
quency dynamics and the a-relaxation peak, which sen- 
sitively shifts with separation from the critical density. 

For a single colloidal particle in a continuum fluid the 
Stokes-Einstein relation connects the particle diffusion 
coefficient and the solvent viscosity, r]D^°^^ = The 
self diffusion coefficient and the mean-squared displace- 
ment at finite colloid densities can, within MCT, be ob- 
tained from the autocorrelation function of the fluctuat- 
ing forces which the single particle experiences from the 
colloidal hquid Q: 

6r^{t) + f di'm^^^^it - t') 5r^{t') = 6D"-^'="t , 

Jo 

(10) 

where ZJ'^ Soif jg ^]-^g short-time diffusion coefficient of 
the single particle which, neglecting HI, is given by 
£)s.Soif = 0. The long-time self diffusion coef- 
ficient, 13^°'^, follows from Eq. ( |lO| ) in the Marko- 
vian limit, D^^-Soif/jjSeif ^ I + dtm^''^^{t) . 

The memory function m^°'^ in MCT is given by an- 
other mode-coupling functional. Thus a priori, within 
MCT one would expect connections or similarities of 
r7i^°'^(w) and G,i{uj) only because of the scaling laws. 
In the /3-scaling region, asymptotically both functions 
exhibit the same shape [|l2|, G'l^{uj)/ha^ x"(w) and 

TO^°'^"(w)//ij„seif x"{^): where x"(^) follows from the 
/3-correlator G(t) in Eq. (^. It is included in Fig. |^. 
The a-superposition principle, Eq. (|^), states that the 
a-relaxation peaks in both functions asymptotically ap- 
proach a density independent shape and shift in paral- 
lel upon varying e. This a-scale coupling also immedi- 
ately predicts the product ZJ^^'^ry to approach a constant 
asymptotically for (j) 4>c fl^- Nevertheless, as for ex- 
ample the a-peak positions need not coincide, the close 
agreement of G^(t) and m^°^^{t) in Fig. |^ over a wide 
window in time or frequency and covering a substantial 
variation in packing fraction is somewhat surprising. Pre- 
sumably it arises, because, during the cooperative struc- 
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tural motion (cage effect), tlie collective density correla- 
tors around the peak in Sq, i.e. on the length scale of 
the average particle distance, dominate the dynamics of 
small-g MCT memory functions. 



D. Comparison with experiments 

The results of the MCT calculations of the previous 
sections, which partially have been testedjn DLS exper- 

can also be 
;3|1 at some- 



JlOl 

18- 



iments aimed at the glass transition |^ 
used to discuss the recent experiments 
what lower densities which were mentioned in the intro- 
duction. Various other aspects of the results and their 
possible experimental relevance have been presented in 
0,|2|,|2|-||, 1^,11] and will not be repeated here. 

As a first aspect, let us point out, that if the mean- 
squared displacement can be measured and thus the con- 
nected memory function, m^'^'^, then the numerical re- 
sults show that a close estimate of the potential part of 
the shear modulus can be obtained. Even beyond 
the connections predicted by the two asymptotic scal- 
ing laws, Eqs. (^,|^), the numerical results exhibited in 
Fig. 1^, show that both functions are closely related, pre- 
sumably because both arise from the cooperative cage 
dynamics. This connection may be considered as a fre- 
quency dependent generalization of the Stokes-Einstein 
relation and was assumed and tested in the recent dif- 
fusive wave spectroscopy measurements of Mason and 
Weitz In another study of the same authors [p3| , 

they also observed that the /^-correlators from Eq. (m 
provide a description of the (directly measured) shear 
moduli spectra in an intermediate frequency window con- 
sistent with the MCT description of the potential part of 

Grj. 

The a-scale coupling predicts that the various relax- 
ation times and transport coefficients of a colloidal sus- 
pension close to the critical packing fraction (p^ shift in 
parallel. For example, the prediction = const, 

for (f) y (j)c follows from Eq. (||). Quantitatively, the 
ratio approaches 5.93 |32j, see Fig. ||, a value very close 
to the classical Stokes-Einstein prediction. Note how- 
ever, that the conditions required for the classical Stokes- 
Einstein relation to hold, clearly are violated at packing 
fractions around the glass transition. A small but no- 
ticeable packing fraction dependence in 1/ {D^'^^^r]) arises 
because of trivial density prefactors connecting the ex- 
hibited moduli of Fig. ^ with the transport coefficients. 
The a-process corrections Eq. (^) and their discussion 
in Eq. (||) suggest that the a-scale coupling should hold 
well for a-relaxation scales obtained at low frequencies. 
This is supported by the observation, that the Stokes- 
Einstein relation considered with limt^^o 4'cG'^{^)/{4'^) 
replacing ry, and nT''^^^"{uj = 0) replacing 1/1?^°'^, consid- 
erably reduces its density dependence. Upon decreasing 
the packing fraction to = 0.9 (/)c (n = 3), where the 
a-scaling law looses validity, this ratio increases by 20 



% relative to the asymptotic value, whereas the actual 
Stokes-Einstein ration increases by 29 %. Even larger 
density dependences can be expected if the long-time dif- 
fusion coefficients are obtained in time or frequency win- 
dows where the dominant corrections to the a-process, 
see Eq. (^, appreciably increase the relaxation times 
of the correlators with shorter asymptotic a-relaxation 
times. This effect is apparent in Fig. |8[ where the wave 
vectors away from the peak position in Sq show an in- 
crease in the relaxation times relative to the asymptotic 
a-scale prediction, which on the other hand holds rather 
well for q = Qp. At q = 1.17qp where Sq{(j)c) — 0-90 a 
65 % increase is seen at n = 3, whereas at q — Qp the 
Stokes-Einstein ratio changes only by 11 %. 

The results concerning the generalized Stokes-Einstein 
relation, which are presented in Fig. H, and their explana- 
tions using Eqs. (^,^ rest on the simphfications caused 
by the bifurcation singularity in the MCT equations and 
by the entailing asymptotic expansions. As discussed, 
the range of validity of the a-scale coupling does not ap- 
preciably exceed e w —0.1 for the studied models of hard 
sphere like colloidal suspensions, and thus microscopic 
corrections to the MCT description of the structural re- 
laxation need to be incorporated in principle beyond this 
distance to (j)c- Thus it is consistent with MCT that for 
packing fractions well below (/)c the coupling of the time 
scales may continue in one system (hard spheres) but not 
in another one (charged spheres) |p4|. Also the estimate 
for the range of packing fractions where a-scaling should 
hold can be expected to be model dependent, as the a- 
master functions, Eq. (||), and their corrections, Eq. (||), 
depend on the fluid structure. 

Approximate expressions might be useful to describe 
the qualitative trends in the intermediate scattering func- 
tions, $g(t), in a wider context, in the same way as the 
Gaussian approximation ||3^ ] is useful for the self inter- 
mediate scattering functions, ^q{t). The Gaussian ap- 
proximation was compared to the MCT results in [ p4[ . 
Recently Segre and Pusey were lead by their DLS scat- 
tering data to propose such a formula and they ob- 
served partial collapse of their rescaled data. 

In Eq. (^) the possibility to collapse the intermediate 
scattering functions onto a common curve is studied us- 
ing the dimensionless time scale Tq '. This is suggested 
by the finding that the a-relaxation amplitudes and re- 
laxation times asymptotically are connected via Tq • see 
Fig. ||. The qualitative connection was expected |1^,|3^] 
but the quantitative closeness surprises and may be pe- 
culiar to the hard sphere system. Satisfactory data col- 
lapse using Eq. (^), see Fig. ^, is possible with deviations 
at long times because of the non-universality of the a- 
process, and at short times, because Eq. (||) violates 
the short -time diffusive motion of colloidal suspensions. 
This reiterates that within MCT there is no connection of 
the obtained long-time diffusion coefficients, which fol- 
low from Eqs. (|l|J^), to the short-time ones The 
only effect of the latter could be a shift in the time scale 
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to , which matches the structural relaxation to the micro- 
scopic transient. 

Similar shapes, however, of the long- and the short- 
time diffusion coefficients with HI were observed by Segre 
and Pusey in the recent DLS experiments on colloid fluids 
below and close to the glass transition ^ . This sim- 
ilarity of the short- and long-time diffusion coefficients 
suggests to collapse the intermediate scattering functions 

with the assumption fl^ : 

exp{-(gV6)p-/i?-Scif),5^2(i)|^ 

which becomes exact for short times. In their experi- 
ments, Segre and Pusey observed data collapse for wave 
vectors starting from somewhat below to the posi- 
tion of the second maximum in Sq. Figure ^ shows the 

solutions of the MCT equations with Z)^ '-^^-' appropri- 
ately rescaled. Reasonable collapse of the curves onto 
a common one, which also is well represented by the 
mean-squared displacement, is observed in a similar wave 
vector range as in the experiments. For short times all 
curves coincide rigorously. For small wave vectors the 
non-diffusive character of the a-process however leads 
to strong deviations for longer times. This trend also 
is present in the experimental data. The non-diffusive 
structural relaxation disagrees with the assumed diffu- 
sive scaling of the density correlators and thus cannot be 
rationalized with considerations of the short-time expan- 
sions following de Gennes. Unavoidable polydispersity 
effects in the experimental data could lead to additional 
deviations from the proposed scaling for small q, but no 
qualitative differences for samples of different polydis- 
persities were reported in ||l^,|2^. Polydispersity effects 
could be incorporated into the present MCT following 
the work for charged colloids in Ref. The partial 

collapse of the data for intermediate and long times re- 
sults from the connections of the a-process amplitudes 
to the time scales discussed in context with the ansatz 
of Eq. (^) . Note that this connection may not be quan- 
titatively satisfied as well in other colloidal systems, like 
e. g. charged colloidal particles Thus the approxi- 
mate scaling may hold less well in other systems. Differ- 
ently from Eq. (j^) , the scaling with the short-time diffu- 
sion coefficients does not rigorously collapse the data at 
a longer time. But the close similarity of t^^ and Dg'-^'-' 
explains that collapse for short times and approximate 
collapse for longer times is achieved. 



IV. CONCLUSIONS 

In the idealized MCT, the long-time dynamics of col- 
loidal liquids is dominated by the structural relaxation. 
Asymptotic expansions close to the critical packing frac- 
tion capture the qualitative aspects of the structural re- 
laxation. In this contribution it is shown that the theo- 
retical results also rationalize some recent experimental 
findings for larger separations from the critical density. 
Corrections to the coupling of the a-relaxation times are 



wave vector dependent as seen in tests of generalized 
Stokes-Einstein ratios . The latter can be generalized 
to finite frequencies as observed in light scattering experi- 
ments , if the potential contribution to the shear mod- 
ulus is considered. The tight coupling of the collective 
density fluctuations as captured in the scaling observed 
by Segre and Pusey on the one hand supports the 
existence of an a-scaling law as predicted by MCT, and 
on the other hand requires the quantitative connection of 
the a-process amplitudes and relaxation times obtained 
here for hard sphere like colloidal particles. The failure of 
the scaling for small wave vectors seen in the experiments 
is predicted by MCT and highlights that the structural 
relaxation cannot be understood from short-time expan- 
sions. Such an approach, often referred to as "de Gennes 
narrowing" concept, would suppose diffusive colloidal dy- 
namics for small wave vectors in disagreement with MCT 
and experiment [ p^ . As the structural relaxation is de- 
termined by the equilibrium structure factor only, hydro- 
dynamic interactions affecting the short-time and tran- 
sient dynamics can be incorporated into the MCT with- 
out changing the long-time predictions. 
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FIG. 1. 

Inverse short-time diffusion coefficients without hydrodynamic interactions (HI), Dq '^'' (left scale), and with HI, 
Dq''^^^ (dot-dashed curve, right scale), versus rescaled wave vector. Without HI, the density variation is determined 
by the structure factor, Do/Dq^^^ — Sq, and is recorded for cf) — (f>c{l — 10^"^^) with n — 1,2,3 (thin dotted, short 
and long dashed line) and (f>c = 0.516 . . . |3^; for n >6, the Sq (almost) collapse onto the bold solid line. Due to the 
rough modeling, the shape of is not varied with density. The inset shows the density dependence of the peak 

position qp for the considered densities corresponding to n — 1,2,3,6,9,12. 

FIG. 2. 

Normalized intermediate scattering functions ^q{t) versus time for five wave vectors indicated by symbols in the 
inset of part (b). In (a) and (b), the height of the circles equals fq/2 and the height of the diamonds gives In 
(a) the packing fraction is (p = 0.999 0c, n ^ 9, and (p = 0.9 0c, n = 3, in (b). Dashed lines marked with circles 
result from calculations without HI and solid lines marked with diamonds from the ones with HI. The inset shows 
the amplitude of the final a-relaxation process. The circles in (a) and (b) indicate the relaxation times estimated 

from <^q{t = Tq) = 

FIG. 3. 

Dimensionless time scales Tq ^ (full circles and left axis) resulting from ff, see Eq. (^), versus wave vector and 
compared to the rescaled a-relaxation times (open circles and right axis) [p6| . The lines through the points indicate 
the corresponding results for the Verlet-Weis Sq from The short dashed curve shows Sq/{q(£f' appropriately 

shifted to match at qp. 



FIG. 4. 

Intermediate scattering functions versus time replotted as suggested by Eq. (^ with Tg^'' taken from Fig. |[ an 
arbitrarily chosen factor enlarges the vertical scale, and the part $q(t) > 0.05 lies in the window for all but one 
correlator. Fig. (a) presents results for a density close to 0c, = 0.999 0c (^^ = 9), whereas Fig. (b) corresponds to 
= 0.9 0c {n — 3). The full solid lines correspond to the mean-squared displacement Sr^lt) from Ref. |3^] scaled 

accordingly. The wave vectors of the exhibited <&q(t) are marked in the inset of part (b). It shows Tq'^ /rj/J (symbols 
and solid line) and the a-relaxation times Tq/Tq^ (dashed line) from Fig. ^. 

FIG. 5. 

Reduced a-relaxation times in units of cf/WODo for various packing fractions and defined by ^q{t = Tq) = fq/2. 
The curves are shifted according to the a-scaling law, Eq. || with 7 = 2.46 [Q, and plotted versus wave vector. The 
thin overlapping lines are for both models at n = 6, 9 and 12. The small circles repeat the rescaled Verlet-Weis result 
from Fig. ^. The bold dashed line results for the model without HI at n = 3 (left scale), and the bold solid line for 
the one with HI( right scale) at n = 3. The inset repeats the data for n = 3 shifted to unity at qp and also includes 
curves for n — 1 and 2 (thin); line styles as in the main part. 



FIG. 6. 

Long-time diffusion coefficients calculated from the times in Fig. |^ and normalized to unity at qp plotted versus 
rescaled wave vector. For n > 6 the thin solid lines overlap and indicate the asymptotic result. Results without HI 
(bold long-dashed line) and with HI (bold solid line) at n = 1 are shown and compared to the rescaled short-time 
diffusion coefficients at n = 1, (thin short dashes) and D^-^™^ (thin dot-dashed line). 



FIG. 7. 

Loss shear modulus, G'/i{uj), solid lines and left scale, compared to the self particle memory function, m^°'^ (lu) 
dashed line and right scale, for the separations from the critical packing fraction corresponding to n = 1, 2, 3, 6 and 
9 as labeled. The inset shows the storage shear modulus for the same packing fractions. The dashed-dotted curves 
in both cases indicate the appropriate Fourier transforms of the /3-correlator from Eq. ^. 
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FIG. 8. 

Packing fraction dependent generalized Stokes-Einstein relations calculated from the viscosity rj and the long-time 
diffusion coefficients from Fig. ^ are shown for various wave vectors indicated in the inset, which shows Sq{4ic)- 
Solid squares use the long-time self diffusion coefficients. The open symbols show the corresponding ratio using 
limtj^o f^cG''(w)/(na;) instead of 77, and m*"^'^ [lo = 0) instead of 1/_D^°'^. 



FIG. 9. 

Intermediate scattering functions with HI effects rescaled using the short-time diffusion coefficients, _Dg from 
Fig. ^ according to Ref. a factor enlarges the vertical scale, and the part ^q{t) > 0.05 lies in the window for 

all but three correlators. Fig. (a) presents Part (a) corresponds to a packing fraction = 0.999 0c {n = 9) close to 
the critical density, whereas part (b) corresponds to a larger separation, (f> — 0.9 0c {n ~ 3). The thick solid lines give 
the mean-squared displacement. The other curves belong to the wave vectors indicated in the inset. It shows the 
normalized times the short-time decay times, r|'^^^^ = l/(q^Dq^^^^) (symbols and solid line) and the a-relaxation 
times Tq/Tq^ (dashed line) from Fig. ^. 



11 




Fuchs and Mayr, Aspects ... Fig. 1 




10"' 10° lO' 10^ 10^ 10' 

tDQm/(f 




10"' 10° lO' 10^ lO"* 10' 



Fuchs and Mayr, Aspects ... Fig. 2 




Fuchs and Mayr, Aspects ... Fig. 3 



I 1 — I I 1 1 llll 1 — I I 1 1 llll 1 — I I I lllll 1 — I I I lllll 1 — I I I lllll 1 — I I I lllll 1 — I I I lllll 1 — I I I lllll 



n = 9 




10 



tDom/d^' 



m — I — 1 1 1 1 iii| — I — 1 1 1 1 iii| — I — 1 1 1 iiii| — I — 1 1 1 iiii| — I — 1 1 1 iiii| — I — 1 1 1 iiii| — I — 1 1 1 iiii| — I — 1 1 1 iiii| 




n = 3 



1.0 - 



0.5 - 



q I Qp 




5 10 qd 



UuX 1 — I I 1 1 llll 1 — I I I I llll 1 — I I I Mill 1 — I I I lllll 1 — I I I lllll 1 — I I I lllll 1 — I I I lllll 



10 10 

Fuchs and Mayr, Aspects ... Fig. 4 
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Fuchs and Mayr, Aspects ... Fig. 5 




Fuchs and Mayr, Aspects ... Fig. 7 




Fuchs and Mayr, Aspects ... Fig. 8 



